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ABSTRACT 

The X-ray emission from normal elliptical galaxies has two major components: soft (kT w 0.2 — 1 
keV) emission from diffuse gas and harder {kT « 6 keV) emission from populations of accreting (low- 
mass) stellar X-ray binaries (LMXB). If the LMXB population is intimately tied to the field stellar 
. population in a galaxy, its aggregate X-ray luminosity is expected to be simply proportional to the 

• optical luminosity of the galaxy. However, recent ASCA and Chandra X-ray observations show that the 

global luminosities of LMXB components in ellipticals exhibit significant scatter (a factor of ~4) at a 

■ given optical luminosity. This scatter may reflect a range of evolutionary stages among X-ray binary 
5^ ' populations in elliptical galaxies of different ages. If so, the ratio of the global LMXB X-ray luminosity 
Qh. to the galactic optical luminosity, Llmxb/ 'L op t, may in principle be used to determine when the bulk of 

stars were formed in individual ellipticals. To test this, we compare variations in Llmxb / L op t for LMXB 
populations in ellipticals to optically-derived estimates of stellar ages in the same galaxies. We find no 
correlation, which suggests that variations in iLMXB/iopt are not a good age indicator for ellipticals. 
Alternatively, LMXBs may be formed primarily in globular clusters (through stellar tidal interactions), 
| rather than comprise a subset of the primordial binary star population in a galactic stellar field. Since 

■ elliptical galaxies exhibit a wide range of globular cluster populations for a given galaxian luminosity, 
' this may induce a dispersion in the LMXB populations of ellipticals with similar optical luminosities. 

[ — ! Indeed, we find that £LMXB/£ op t ratios for LMXB populations are strongly correlated with the specific 

t— I ■ globular cluster frequencies in elliptical galaxies. This suggests that most LMXBs were formed in globular 

clusters. If so, Chandra observations of central dominant galaxies with unusually large globular cluster 

populations should find proportionally excessive numbers of LMXBs. 

, Subject headings: binaries: close — galaxies: elliptical and lenticular — globular clusters: general — 

■ X-rays: binaries — X-rays: galaxies 

■ 

i ' 1. introduction fraction of stars in binaries, the distributions of binary sep- 

mi v c i it i i i arations and mass ratios, etc.) are largely independent of 

$H I he X-ray emission from normal elliptical galaxies has ; , . . , . , , , , , r T * ri-T-. n ,, . 

-4— > ■ . c, /, m n o i i ?n • • their galactic context, the number ol LMXBs hence their 

r/3 , two maior components: soft (kl w 0.2 — 1 keV) emission , v i • \ ■> 111 • i , • ^ 

d ■ from diffuse gas and harder (kT « 6 keV) emission from aggregate X-ray luminosity) should be simply proportional 

• • , , ,. c . • n \ 4- ii v u- to the number of stars m the galaxy (and thus their total 

^ populations of accreting (low-mass) stellar X-ray binaries ... . . , ° J y 

• i-h (LMXB). The X-ray properties of the LMXB component P TT . , . , ,. ™ > i 

k> n , j-rc i, , j , j ,i ■ i r High angular resolution Chandra observations are now 

JC, have been difficult to determine, due to their spatial coniu- ., . ..... . , , . . ril ,. r 

r 5_j . .,, ,.rn . . , , i it j allowing individual LMXBs to be resolved out of the dit- 

« sion with diffuse gaseous emission and to spectral nardness r ° „ ... , ... . 

which places much of the LMXB emission outside the effec- ^ S ?000^1 -1 " IrwTn TbI man 2000 S 

tive bandpasses of most imaging X-ray satellite spectrom- e . a ',. . ' ' ; arazm, rwin, regman , , 

, rrnl r j t tmvd i i l Angclmi, Locwcnstcm, & Mushotzky 2001: rmogucnov & 

eters. Ihe presence ol the LMXB component has been T b „„', ... ' ,. . J . ' , ° . . 

■ r j- j ji tt c 4-itj Jones 2001), which makes their composite spectral analy- 

mferred m part through observations ol spectral harden- . . , „ ,. . -v. ... , 

„. ,. , -.t -l n v i 4-- i sis much easier. Ihe bulk ol the hard emission m normal 

mg m ellipticals with progressively smaller X-ray to optical r T * r ^-r^ i i r 

, . ., ,. t-,,. f, , , . n m • i • ■ lnnn i • ellipticals indeed comes from LMXBs, rather than from 
luminosity ratios (Kim, labbiano, & Irmchien 1992), m- . 1 . . . ; . n ' . , 

i • i • 1 1 i 1 1 i , , ,.,, n . ,i advcction-dommatcd accretion hows onto massive, central 

dicatmg that they have relatively little gas, exposing the . . . , , , , ... ,. _ , „ _ , . ' 

i i t a /tvd -t n i a- uiuvD i black holes, proposed by Allen, di Matteo & rabian (2000 . 

harder LMXB component. Populations ol LMXBs are also TT , , J , ,,. . . , , 

, j • it i i i i ■ i v j j i ■ Until Chandra observes more nearby ellipticals, the 

expected m ellipticals simply by analogy with detections . . , , . , , . , . ' . 

r j. , . v. u -j i a strongest spectral constraints to date on the hard stellar 

ol discrete LMXB sources m nearby spheroids such as the T1 . T „° r 

, , j i i r /-ii j atoi ("n t LMXB component m a large sample oi ellipticals will still 

bulges and halos ol our Galaxy and M31 (lorman, Jones . ,„„. , , , , , T Tn •, 

p m. i moc n • P .i. p mL- u- ■ mo7\ come from AtiCA spectra (Matsumoto et al. 1997; White 

& lucker 1985: Camzares, rabbiano & Irmchien 1987), ™™ „„„, r,n^r,N ,, rr ,• i , p 

as well as in the radio tralaxv Ceutaurns A TTnruer et a 2000 ' 2001 ' 2002 )' SmCe the effect ive angular resolution of 

as well as in Liic i aciio yaiaxy v^ciiiaui us i _luiiici ei ai. . ... „ nl lr .. 

,qq 7 N AtiCA imaging spectrometers is 2 — 3 (hall-power diam- 

a . i , . ,i . ;t . , i v i eter), confusion prevents individual LMXBs from being 

A simple argument suggests that the total X-ray hum- . . , „ . ™ 

•4-- tTHVD i j? • it a- i ■ v.i t easily resolved out oi the diffuse gas m ellipticals. Ihe 

nosities ol LMXB populations m ellipticals might be pro- . , J T%rvr> , , & , „ ,. ,. . , , 

,. i. n i ,,, , . •/,, hard LMXB component can be spectrally distinguished 

portional to the stellar luminosities ol the galaxies: it the . n v , . , 

r , , n . , , Irom the sotter gaseous component, however. Matsumoto 

properties oi low mass binary stellar systems (such as the ° 
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et al. (1997) separated the hard LMXB component in ellip- 
ticals from softer gaseous emission by considering ASCA 
GIS spectral energies above 4 kcV. Stacking 4-10 keV GIS 
spectra for 12 ellipticals, they found a best-fit thermal 
bremsstrahlung model with kT = 12^ 9 ^ keV (where er- 
rors arc 90% confidence limits) ; a power-law with photon 
index a = 1.8+o' 4 fit equally well. Matsumoto et al. (1997) 
found that the X-ray luminosities of the LMXB compo- 
nents were proportional to the optical luminosities of the 
ellipticals, but with a surprisingly large scatter (a factor 
of ~ 4). Some of the ellipticals included in the sample 
have significant X-ray emission from active galactic nuclei 
(AGN), which may account for some of the scatter. 

White (2000, 2001, 2002) performed a similar ASCA 
analysis, but with a larger spectral bandwidth (0.7-10 
keV), on six normal ellipticals (i.e., ellipticals without 
significant AGN emission). Spectra were extracted from 
within metric radii of 6r e from the galactic centers, where 
r e is the optical effective radius of a galaxy. The GIS 
spectra of the six ellipticals were simultaneously fit with 
both soft (gaseous) and hard (LMXB) emission compo- 
nents. (Only GIS data were used because the GIS de- 
tectors have twice the effective area of the SIS detectors 
above 7 keV.) The temperatures (or power-law indices) of 
the hard components in the galaxies were tied together, 
while the temperatures of the soft components (if present) 
were allowed to vary individually. Much tighter spectral 
model constraints were provided by the increased spectral 
bandwidth compared to the 4-10 keV bandwidth in the 
stacked spectral study of Matsumoto et al. (1997). The 
spectra of the LMXB components were fit equally well 
by a bremsstrahlung model with kT — 6.4^ ^ keV or a 
power-law model with photon index a = 1.82+°;^ (errors 
are 90% confidence limits). Individual fits to each galaxy 
in the set were consistent with the results of the joint fits 
and fluxes were obtained by adopting the best jointly fit 
temperature. These are the tightest constraints to date on 
the global spectral properties of the stellar LMXB compo- 
nent in ellipticals. GIS X-ray fluxes were determined for 
the LMXB components in an additional 8 ellipticals which 
had poorer photon statistics by fixing the temperature of 
the hard component to 6.4 keV. The resulting X-ray lumi- 
nosities of the LMXB components in the 14 galaxies were 
found to be proportional to the optical luminosities of the 
galaxies, with a factor of ~ 4 — 10 scatter. 

Recent Chandra observations (Sarazin et al. 2000, 2001: 
Angelini et al. 2001; Blanton, Sarazin, & Irwin 2001) 
also show that the X-ray luminosities of the resolved 
LMXB components in ellipticals exhibit significant scat- 
ter at a given optical luminosity. (The composite spectra 
of these LMXB components are also consistent with the 
hard ASCA spectral results described above.) Although 
their scatter in luminosity is much smaller than that of 
the softer gaseous component (which ranges a factor of 
100 in X-ray/optical flux ratio), it is still larger than ex- 
pected if the LMXB component is strictly proportional to 
the stellar component. Following a preliminary study by 
White (2000, 2001), we next consider two possible sources 
of variance in the X-ray luminosities of LMXB populations 
in elliptical galaxies: age differences among ellipticals and 
differences in their globular cluster populations. 



2. POSSIBLE SOURCES OF VARIANCE IN LMXB 
POPULATIONS IN ELLIPTICALS 

2.1. Variations in Galaxy Ages 

If LMXB populations are a subset of the primordial bi- 
nary systems in a galactic stellar field, then the evolution 
of an LMXB population should be tied to the evolution 
of visible stars. The evolution of the aggregate X-ray lu- 
minosity of an LMXB population is best assessed with a 
population synthesis calculation (Wu 2001). A simpler, 
more phenomenological approach was adopted by White 
& Ghosh (1998). The luminosity of an LMXB population 
(neglecting burst sources) is expected to rise for about a 
Gyr (White & Ghosh 1998), driven by the time it takes the 
less massive secondaries in binaries to evolve off the main 
sequence, overflow their Roche lobes, and start dumping 
mass onto the more massive (primary) compact stellar 
remnants in the binaries. After this Gyr ramp-up, the 
luminosity of an LMXB population is expected to slowly 
decline (White & Ghosh 1998; Wu 2001). 
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Fig. 1. — Mean ages of the stars in elliptical galaxies (Trager et 
al. 2000) are plotted against the ratios of the X-ray luminosity of 
their LMXBs to their optical B-band luminosity (White 2002). The 
open circle indicates an upper limit to Llmxb due to the presence 
of an AGN. 

The optical stellar luminosities of galaxies are also 
slowly declining with time. Figure 1 of Bruzual & Char- 
lot (1993) shows that the (visual) optical luminosity of 
an elliptical galaxy declines as Ly cx i~ - 95 . This time 
dependence is close enough to estimates for the temporal 
decline in the X-ray luminosity of LMXB populations that 
it is difficult to assess how the X-ray to optical ratio of the 
LMXB component is evolving at the present epoch. 

For our purposes, ambiguity in the theoretical expecta- 
tion for the temporal evolution of iLMXB/^opt does not 
matter — we can deduce the evolution empirically. If vari- 
ations in the X-ray properties of the LMXB component 
for galaxies of a given optical luminosity are correlated 
with independent age estimates, the Llmxb/ Lopt ratio it- 
self might then be used as an indicator of galactic age. 

To test this possibility, we compare the £ L mxb / L op t ra- 
tio for the LMXB components from the elliptical sample 
of White (2002) to recent stellar age estimates for ellip- 
tical galaxies (Trager et al. 2000). The stellar age esti- 
mates were derived by comparing the strengths of opti- 
cal stellar spectral line indices of <Fe>, Mg b and H/3 
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(which were corrected for the effects of non-solar abun- 
dance ratios), to the single- burst stellar population mod- 
els of Worthcy (1994). For their sample of 50 ellipticals, 
Trager et al. (2000) find ages ranging from 1.5 — 18 Gyr. In 
Figure 1, the relevant optical stellar age estimates (rang- 
ing from ^4 — 14 Gyr) are compared to the Llmxb /-^opt 
ratios for the LMXB populations in the eight ellipticals of 
(the fourteen in) White (2002) that were included in the 
Trager et al. (2000) sample. The X-ray luminosities were 
calculated in the 0.5-4.5 keV range and their errors are 
90% confidence limits. There is apparently no correlation, 
so variations in the Llmxb /L opt ratios for LMXB popula- 
tions are not likely to be due to galactic age differences. 

2.2. Globular Cluster Population Variations 

Alternatively, most LMXBs may be formed in globular 
clusters, rather than in the stellar field of a galaxy. Obser- 
vations of LMXBs in our own Galaxy show that it is very 
difficult to make LMXBs in the field. There arc compa- 
rable numbers of high mass X-ray binaries (HMXBs) and 
LMXBs, but HMXB lifetimes are only ~ 10 5 yr, which 
are 10 2-5 times shorter than LMXB lifetimes. Thus, the 
formation of LMXBs in our Galaxy is 10 2 ~ 5 times less 
efficient than the formation of HMXBs. 

Meanwhile, globular clusters are the most efficient for- 
mation site for LMXBs in our Galaxy: ~ 20% of all 
LMXBs in our Galaxy reside in globular clusters, yet glob- 
ular clusters contain < 0.1% of the stars in our Galaxy 
(Katz 1975). Thus, globular clusters form LMXBs > 200 
times more efficiently than the rest of the Galaxy, through 
tidal capture during close stellar encounters (Clark 1975). 

Given this, it is reasonable to consider whether most 
LMXBs are formed in globular clusters. Grindlay (1984) 
suggested that X-ray bursting binaries in the bulge of our 
Galaxy (not residing in globular clusters) were formed in 
globular clusters which were later destroyed by tidal shocks 
while passing through the Galaxy's disk. This disruption 
mechanism is not relevant to globulars residing in ellipti- 
cal galaxies, but if globulars are on eccentric orbits, they 
may still lose stars due to galactic tides at perigalacticon. 
Vesperini (2000) estimates that only ~ 10% of globular 
clusters are disrupted in the most luminous ellipticals, but 
the fraction disrupted increases considerably for progres- 
sively smaller elliptical luminosities. Also, some LMXBs 
formed in globular clusters may be ejected by the velocity 
kicks imparted by supernovae during the formation of neu- 
tron stars. Typical supernovae kick velocities are thought 
to range from 100 — 500 km/s (Kalogera 1996; Terman, 
Taam, & Savage 1996; Fryer, Burrows, & Benz 1998), 
while the escape velocities of globular clusters range from 
5-60 km/s (Webbink 1985). Since LMXB lifetimes can be 
quite long (> 10 9 yr), LMXBs ejected from globular clus- 
ters can be visible outside the clusters for a long time. If 
LMXBs are made primarily in globular clusters, we would 
expect their number to be proportional to the number of 
globular clusters, regardless of the fraction ejected from 
clusters. 

With the excellent angular resolution of Chandra, it 
is possible to determine positions of LMXBs with suffi- 
cient accuracy to determine if they coincide with globular 
clusters. Observations of ellipticals show that a signifi- 
cant fraction are located in globulars (Sarazin et al. 2000, 



2001; Angclini et al. 2001; Randall, Sarazin, & Irwin 2001; 
White, Davis & Hancs 2002). The fraction ranges from 
>20 to 70%. These are lower limits, since deeper Chandra 
exposures may detect fainter X-ray sources in additional 
globulars. In the same region of these galaxies, the glob- 
ular clusters typically provide only <~ 0.1% of the optical 
light, which implies that optical stars in globular clusters 
are much more likely to be the donor stars in LMXBs 
than in the field. The primary limitation with identifying 
LMXBs with globulars in ellipticals at present is the lack 
of Hubble Space Telescope (HST) globular cluster lists for 
many of the Chandra observed ellipticals. 

Since elliptical galaxies exhibit an order of magnitude 
range of globular cluster populations for a given galaxian 
luminosity, we can assess whether these variations are cor- 
related with variations in L^MXB/L op t- In Figure 2, we 
compare the specific frequency of globular clusters S (the 
number of globular clusters per galaxy visual luminosity 
in units of the luminosity of a My = —15 galaxy) from 
the compilation of Kissler-Patig (1997) to the iLMXB/^opt 
ratio for the LMXB populations in the elliptical sample 
described above. The Kissler-Patig (1997) compilation 
includes error estimates for S, as indicated in the fig- 
ure, for all but one elliptical in this sample. There ap- 
pears to be a strong correlation (White 2000, 2001), and 
the relationship is consistent with direct proportionality: 
LhMXB/Lopt x S 11 - 2 * 4 . Thus, LMXB populations may 
indeed be controlled by globular cluster populations. 
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Fig. 2. — Logarithm of the specific frequency of globular clusters S 
(Kissler-Patig 1997) vs. the X-ray/optical luminosity ratio for ellip- 
tical galaxies. The best fit relationship (using the ASURV statistical 
package (Isobe et al 1986) is Llmxb /£ op t oc S 1 - 2 * ' 4 , with a 90% 
probability of correlation via Kendall's T test. The open circle in- 
dicates an upper limit to Llmxb due to the presence of an AGN. 



3. DISCUSSION AND CONCLUSIONS 

This study suggests that globular cluster populations 
in ellipticals control the evolution of their populations of 
low mass stellar X-ray binaries. At the present epoch, age 
differences among the ellipticals in the small sample de- 
scribed above do not seem to be correlated with the varia- 
tions in their X-ray binary populations. It would be useful 
to compare the LMXBs in early-type galaxies with more 
detailed properties (such as color and magnitude) of their 
globulars (other than just the number), to see if the history 



4 



WHITE, SARAZIN, & KULKARNI 



of star formation in the globular affects the LMXB popu- 
lation. The Chandra observations of NGC 1399 suggested 
that LMXBs are more likely to be found in redder globular 
clusters (Angelini et al. 2001). There also is evidence that 
LMXBs are more likely to be found in brighter globulars 
(Angelini et al. 2001), although it is unclear whether this 
is just due to the larger number of stars in brighter glob- 
ulars, or indicates a higher probability per star in more 
massive clusters. 

It is possible that most or all LMXBs in early-type 
galaxies are formed in globular clusters. In this scenario, 
any field LMXBs in ellipticals would have escaped from 
globulars. The field LMXBs might have been ejected by 
the kick velocities resulting from supernovae, by stellar dy- 
namical processes such as superelastic encounters in which 
internal binding energy in the binaries is converted into ki- 
netic energy of motion of the binary, or by the dissolution 
of the globular due to tidal effects. 

There is some evidence that the LMXBs in globular clus- 
ters in elliptical galaxies are brighter in X-rays than those 
in the field (Angelini et al. 2001). If all of the LMXBs 
were formed in globular clusters, this would require that 
the less luminous and presumably less massive LMXBs 
systems were preferentially ejected from globular clusters. 
Alternatively, if a significant fraction of the field LMXBs 
were formed in the field, this would imply differing stellar 
and binary evolution histories in the field and in globular 
clusters. 

Many sources associated with globulars have X-ray lu- 
minosities which exceed the Eddington luminosity of a 1.4 
Mq neutron star, LEdd,NS ~ 2 x 10 38 ergs s _1 (Sarazin et 
al. 2000, 2001; Angelini et al. 2001; Blanton et al. 2001; 
Randall et al. 2001; White et al. 2002). In general, LMXB 
luminosity functions within early-type galaxies have a 
break at a luminosity very close to LEdd,NS (Sarazin et al. 
2000, 2001; Blanton et al. 2001; Randall et al. 2001). This 
break may separate accreting black holes, at higher lumi- 
nosities, from accreting neutron stars at lower luminosities 
(Sarazin et al. 2000). If the sources are Eddington- limited 
single binaries, then the brightest sources associated with 
globular clusters must contain quite massive (> 20M Q ) 
black holes (Sarazin et al. 2000). It is not clear how globu- 
lar clusters would form and retain such massive black holes 
in LMXBs and also not have them form binary black hole 
systems (e.g., Portegies Zwart & McMillan 2000). 

On the other hand, these super-Eddington sources might 
be globular clusters containing multiple X-ray sources (An- 
gelini et al. 2001). The difficulty with this suggestion is 
that the fraction of globulars which contain at least one X- 
ray source is low; the fraction is about 3-4% in NGC 1399 
(Angelini et al. 2001), NGC 4697 (Sarazin et al. 2000, 
2001), and in NGC 4649. (It is interesting that this frac- 
tion is almost the same in these three ellipticals.) The 
luminosity function of the LMXBs isn't steep enough that 



multiple sources should affect the luminosity function at 
the high end, given this low fraction. On the other hand, 
if the LMXBs are mainly confined to a smaller subset of 
globular clusters (e.g., very bright and/or very red and/or 
core-collapsed systems), then the fraction of the relevant 
subset of globulars with X-ray sources and the probability 
of multiple sources might be much higher. 

There is some evidence that these luminous super- 
Eddington X-ray sources are associated with globular clus- 
ters in elliptical and SO galaxies, but not with the bulges of 
spiral galaxies. There are no super-Eddington LMXBs as- 
sociated with globular clusters in the bulge of our Galaxy 
(e.g., Hut et al. 1992). The same is true of the central 5' 
of the bulge of M31 (Shirey et al. 2001). The luminosity 
of LMXBs in the nearby large spiral bulge NGC 1291 has 
a cut-off at a luminosity which is very close to LEdd,NS- 
Unless this is the result of statistical fluctuations and the 
large populations of globular clusters in elliptical galaxies, 
this suggests that the binaries in globulars in ellipticals 
have different histories than those in spiral bulges. 

In general, the specific frequency of globular clusters 
S increases along the Hubble sequence from late to early- 
type spirals, from spirals to SOs, from SOs to ellipticals, and 
from normal giant ellipticals to cDs (e.g., Harris 1991). If 
LMXBs all are born in globular clusters, then one would 
expect the specific frequency of LMXBs and their X-ray to 
optical ratio to increase along the same sequence. In par- 
ticular, a number of central dominant ellipticals in galaxy 
clusters are notable for having enormous populations of 
globular clusters. White (1987) suggested that these un- 
usually populous systems belong to the galaxy clusters, 
rather than individual galaxies: globular clusters will have 
been tidally stripped from individual galaxies during the 
collapse of a galaxy cluster and will then virialize along 
with other cluster constituents. If LMXBs are formed pri- 
marily in globular clusters, then central dominant galax- 
ies with high specific globular cluster frequencies (such as 
M87, NGC 1399, NGC 3311 and NGC 4874) should be 
particularly rich in LMXBs, as well. Chandra observations 
have already confirmed the large populations of LMXBs in 
NGC 1399 (Angelini et al. 2001), and in the region of this 
galaxy with a globular cluster list based on HST observa- 
tions, 70% of the LMXBs are located in globular clusters. 
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